Excitons and Optical Absorption on the Surface of a Strong Topological Insulator 

with a Magnetic Energy Gap 
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We present a theoretical study of interacting electron-hole pairs located on a magnetized surface 
of a strong topological insulator (TI). The excitonic energy levels and the optical absorption on 
such surface display unique and potentially measurable features such as (i) an enhanced binding 
energy for excitons whose total angular momentum is aligned with the magnetic exchange field, (ii) 
a stark dependence of the optical absorption on the direction of the magnetic exchange field as well 
as on the chirality of the incident light, and (iii) a tunable center-of-mass motion of spinful excitons 
induced by particle-hole asymmetry in the exchange field. Our predictions are relevant to surfaces of 
magnetically doped TIs or surfaces coated with magnetic films, in addition to TI nanowires placed 
under longitudinal magnetic fields. 

PACS numbers; 



I. INTRODUCTION 

The discovery of strong topological insulators has cul- 
minated in a new classification scheme for solids»i Strong 
TIs are bulk insulators endowed with a topological invari- 
ant that manifests itself through robust metallic surface 
states containing a nontrivial Berry phase. These sur- 
face states have inspired numerous attractive proposals 
for spintronics applications^ because they are strongly 
spin-orbit coupled and because they exhibit essential ro- 
bustcness against non-magnetic perturbations. In addi- 
tion, they have a linear energy-momentum dispersion at 
low energies, they may display a half-quantized quantum 
Hall effect when time-reversal symmetry is broken, and 
they may host exotic particles such as Majorana fcrmions 
when placed in proximity to a superconductor. 

Thus far both experimental and theoretical studies of 
topological surface states have either ignored electron- 
electron interactions or treated them ad hoc through 
mean field^ theories. This omission is partly justified 
by the fact that weak, short-range interactions between 
massless Dirac fcrmions arc irrelevant in the renormal- 
ization group (RG) sense. Experimental support for the 
RG argument originates from the sizeable Fermi veloci- 
ties {vf — 5 X lO^m/s) and the large static dielectric con- 
stants (e ~ 50—100) found in several well-established TIs, 
which considerably reduce the value of the bare Coulomb 
coupling g = [e^ / e) / {2hv p) {g ^ 0.05). However, there 
appears to be no fundamental principle that forbids the 
existence of TIs with smaller— dielectric constants and 
thereby more strongly interacting surface states. This 
expectation is made more plausible by the ongoing and 
rapid expansion of experimentally confirmed TIs, which 
currently encompass a variety of bismuth compounds^ as 
well as chalcogenides^. For instance TlBiSe2, recently^ 
identified as a TI, has been previously reported^ to dis- 
play a relatively small static dielectric constant of e ~ 25 
in thin films. Likewise HgTe, which has been recentlj^ 
coaxed into a strong topological insulating state, has 



e ~ 20 in bulk. Moreover, the irrelevance under RG 
of electron-electron interactions is questionable when the 
chemical potential of the topological surface states is lo- 
cated inside a magnetically generated energy gap (which 
introduces an infrared cutoff in the RG fiow equations). 

This paper aims to initiate a systematic study of the 
interplay between spin-orbit interactions, broken time re- 
versal symmetry, and Coulomb interactions on surfaces 
of strong TIs. In particular, we focus on low-temperature 
optical properties of topological surface states whose 
chemical potential lies inside a magnetically generated 
energy gap. The underlying microscopic mechanism for 
this energy gap at the surface is an exchange field in 
the direction perpendicular to the surface, which can be 
produced by magnetic doping^ or else by deposition of 
a ferromagnetic insulator with perpendicular anisotropy 
(the latter is yet to be realized experimentally). Of 
course, a magnetic gap may also be opened by the ap- 
plication of a perpendicular magnetic field; however, in 
this case the surface states are modified qualitatively ow- 
ing to the formation of Landau levels. Herein we disre- 
gard orbital effects^ and thus limit ourselves to exchange- 
induced magnetic gaps. 

The simplest interacting problem one can conceive in 
this playground is that of a surface electron-hole excita- 
tion across the magnetic bandgap. Some of the peculiar- 
ities of this two-body problem are made apparent in Sec- 
tion II. On one hand, the relative motion of the electron- 
hole pair is inextricably linked to its center-of-mass mo- 
tion, which is conducive to phenomena like the spin-Hall 
effect of excitons. On the other hand, optically generated 
electron-hole pairs have a nonzero center-of-mass group 
velocity provided that the exchange field is particle-hole 
asymmetric. While the former property is also present 
in certain ordinary semiconductors)^ the latter trait is 
exclusive of topological surface states. 

In Section III we consider a two dimensional surface 
with translational and rotational symmetry. In this case, 
broken time reversal symmetry and spin-orbit coupling 
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conspire to produce Coulomb interaction matrix elements 
that favor a nonzero total angular momentum for the 
electron-hole pairs. This skeweness of the Coulomb in- 
teraction is transferred to the binding energy of dilute 
excitons, which becomes largest for the p-wave angular 
momentum channel. When the surface is irradiated with 
linearly polarized light, this p-wave exciton is optically 
bright and appears prominently in the absorption spec- 
trum of systems with g > 0.05. In contrast, for circularly 
polarized light the dominant p-wave exciton is optically 
dark or bright depending on the handedness of the light 
as well as on the orientation of the exchange field. 

In Section IV we turn to low dimensional TIs. Typ- 
ically nanowires show more pronounced excitonic and 
electric-field effects in their optical absorption spectrum. 
We verify that, in presence of particle-hole asymmetry, 
optically generated excitons have a nonzero center-of- 
mass group velocity. Incidentally, we establish the ap- 
pearance of electric-field-induced Franz-Keldysh oscilla- 
tions. Nanostructures made of TIs are being routinely 
fabricatedfii which may motivate optical measurements 
that test the effects discussed herein. 

Our quantitative results contain unknown parameters 
such as the surface dielectric constant and the lifetime 
of the excitons; it may be possible to extract these 
from photoluminescence experiments that are presently 
under wayi^ in TIs. 

II. GENERAL CONSIDERATIONS 

In this section we consider an interacting electron-hole 
pair located at the surface of a magnetized STI, and dis- 
cuss some of its generic aspects. The Hamiltonian de- 
scribing this problem is 

H = VpCTe ■ {z X Pe) + AgfJ • (Tg 

- vpCTh ■ [z X ph) + A,,ri -ah ef • r, (1) 

er 

where is the Fermi velocity, z is the unit vector nor- 
mal to the surface, ^e{h) — is the momentum 
of the electron (hole), re(/j) is the position of the electron 
(hole), r = Tg — Tfi is the relative coordinate, cTg = cr ® 1 
is the spin operator acting on the electron (<t is a vec- 
tor of Pauli matrices and 1 is the 2x2 identity matrix) , 
and <T/i = 1 (X) cr is the spin operator acting on the hole. 
Ae(/j) = Ae(/j)f2 is the exchange field sensed by the elec- 
tron (hole), which can be generated either by magneti- 
cally ordered impurities doped in the STI bulk, sprinkled 
on the surface^ or else by deposition of a ferromagnetic 
insulator overlay. D, denotes the direction of the mag- 
netization and is assumed to be spatially uniform. The 
exchange field originates from virtual transitions of the 
surface fermions onto localized magnetic ions and need 
not be the same for electrons and holes (i.e. A/;). 
The reasons for this asymmetry are twofold. On one 
hand, the exchange (or Kondo) coupling between itiner- 
ant fermions and localized magnetic moments is gener- 



ally energy- dependentr^ i.e. A = A(e) where e = hvpk. 
On the other hand, the magnetic impurities are gener- 
ally not particle-hole symmetric, i.e. adding an electron 
to the local moment costs different energy from subtract- 
ing one. It follows that hm^^o dA/de ^ 0. In Eq. ([T]) 
we adopted A(e) = Aee(e) + AhG{-e) (where 6(2;) is 
the step function), which constitutes the simplest real- 
ization of broken particle-hole symmetry in the exchange 
coupling. While often ignored for single-particle physical 
observables, Ag 7^ A^, may lead to potentially observ- 
able two-body effects as indicated below. Throughout 
this work we assume il^ 7^ 0, which creates an energy 
gap of (Ae -|- A}i)^z for non-interacting particle-hole ex- 
citations. In addition, we assume that the Fermi energy 
of the surface states is located within this energy gap. 
The attractive Coulomb interaction between the electron 
and the hole is screened by the static dielectric constant 
e, whose value is expected to be a weighted average of 
the dielectric constants corresponding to the bulk TI and 
the vacuum (or the ferromagnetic insulator). The larger 
weight is ascribed to the medium into which the surface 
states decay more slowly. £ is the in-plane electric field 
acting on the surface states. 

Adopting center-of-mass (COM) and relative coordi- 
nates, Eq. (m can be rewritten as 

n = VF {(Te - (Th) • X P + (2) 

+ {f^e + (Th) ■ ( z X p + —Cl] e£ -r, 

2 \ fi? / er 

where p = Pe — p/i is the relative momentum, P = (pe + 
Pft)/2 is the COM momentum and A± = (Ae ± A,i)/2. 
Ce + (Th is the total spin operator for the electron-hole 
pair, and it couples to the relative momentum. This can 
be easily understood by recalling that, for a single Dirac 
fermion, the spin of an electron in the "conduction" band 
with momentum k is perfectly antialigned with the spin 
of an electron with the same momentum in the "valence" 
band. The in-plane component of the exchange field can 
be absorbed by a shift in the relative momentum of the 
electron-hole pair, and therefore it is inconsequential. 

Because both electric fields and Coulomb interactions 
couple only to the relative coordinate, the COM momen- 
tum P is a good quantum number. Within the enve- 
lope function approximation, the conservation of P en- 
ables a factorization of the electron- hole wavefunction via 
*(R,r) = exp(iP • R)x(P,r), where R is the COM co- 
ordinate and X is the "relative" wave function. The co- 
dependence of X on P and r hints at the non-separability 
between the center-of-mass motion and the relative mo- 
tion. In effect, the COM and relative parts in Eq. ([5]) 
(i.e. the first and second lines) do not commute with 
one another. This observation applies to graphene and 
carbon nanotubes as wellf^ however, unlike in graphene 
and carbon nanotubes, in TIs momentum is coupled with 
spin rather than with a (less manipulable) sublattice de- 
gree of freedom. The coupling between COM and relative 
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degrees of freedom creates prospects for spintronics ap- 
plications, such as a spin-Hall effect of excitonsiS derived 
from the Berry connection Aj ~ i J drx*dx/dPj. 

Aside from P, the 2-component of the electron-hole 
pair's total angular momentum is also a good quantum 
number so long as P = £ ^ A_ = 0. Indeed one can 
readily verify that limp^g ^a_-!.o T^^^'HT?. = "H for 7?. = 
exp(is^0) exp(iP0), where — + cr^, P — xpy — ypx 
and (j) is an arbitrary angle of rotation. In Section III we 
shall take advantage of this U(l) symmetry in order to 
simplify the calculation of the excitation spectrum. 

An intriguing peculiarity of topological surface states is 
that particle-hole asymmetry in the exchange field pro- 
duces a fictitious electromagnetic vector potential a = 
/S.-/vfz X VI for the center-of-mass coordinate. There- 
fore, the mechanical or physical COM momentum of the 
electron- hole pair, 11 = P — a, is different from the canon- 
ical momentum P = —ihd-R,. Accordingly, optically ex- 
cited electron-hole pairs (which are essentially vertical 
across the bandgap and therefore have P = 0) contain 
a nonzero mechanical momentum. As we show in Sec- 
tion IV, the group velocity of these excitons is given by 
Vg ~ vf{^- /^z^+)z X Cl for vf\sl\ <C r2zA+ (which 
implies Vg << vp), and consequently their motion may 
be controlled by tilting the direction of magnetization 
away from z. Admittedly A_ itself is not easily tunable 
in a single topological surface. However, nanoribbons 
with two Dirac cones (one on each surface) offer a plat- 
form wherein A_ is tunable (by changing the relative 
magnetization of the two surfaces) for spatially indirect 
electron- hole pairs. In this paper we limit ourselves to 
topological surfaces with a single Dirac cone. 

III. ELECTRON-HOLE PAIRS ON A 
TWO-DIMENSIONAL SURFACE 

This section is devoted to the characterization of opti- 
cally generated excitons and their photolumincsccnce on 
a two dimensional topological surface with a magnetic 
energy gap. For simplicity we assume Ae = Ah = A 
and £ = 0. Moreover we take A = Az, which involves 
no loss of generality because the in-plane component of 
the exchange field can be shifted away by redefining a 
new origin in momentum space. These choices render the 
problem more easily soluble in momentum space, where 
the low energy effective Hamiltonian may be written as 
n='Ho+ll, with 

a—c.v k 

a k,k',q 

+ 12 E ^^""'^U+q^U'-q^-^k'^Ck, (3) 
k,k',q 



i?k,c = —Ek,v = Ek = \/h'^vjpk^ + A^ is the energy dis- 
persion for non-interacting massive Dirac fcrmions (mass 
cx A), and 

Fj^' = (a,k + q;a,k'-q|Z^|a,k;a,k') 

F'^' = (c,k + q;w,k'-q|Z^|c,k;i;,k') (4) 

are the intraband and interband Coulomb matrix el- 
ements. The operators that create Bloch eigenstates 
in the "conduction" and "valence" bands (hereafter de- 
noted as |ck) and jt'k)) are ^'^ ^ = cos{9\^/2)cl^^ + 
exp(z(pk) sin(6'k/2)c[^ and = sin(6'k/2)4_.,. - 

cxp(i</7k) cos(0k/2)cjj. ^, respectively, ^ creates a free 
(in vacuum) electron with momentum k and spin a. Also, 
cpk = tan~^ (ky/kx) +tt/2, cos6'k = A/Ek and is the 
area of the topological surface. 

Let us consider a single clcctron-holc pair, where the 
electron is in the conduction band and the hole is in the 
valence band. Taking P = 0, the wavefunction for this 
electron-hole pair is 



|xc)=^x(k)<k*-,k|0), (5) 

k 



where |0) is the ground state without particle- hole exci- 
tations, namely completely filled valence band and com- 
pletely empty conduction band. x(k) can be found 
from the requirement that |xc) be an eigenstate of H. 
?^|xc) = eph|xc) is satisfied if and only if 



{2Ek + <5Efe)x(k) -^T. ^C'' l^')^(k') = ephx(k), (6) 

k' 



where V(k,k') = (c, k; u, k'|W|c, k'; u, k) (27reVe|k - 
k'|)(c,k|c,k')(w,k>,k). In the derivation of Eq. dH) 
we have treated the interband Coulomb interaction 
(]/(^)) exactly while adopting a mean-field Hartree-Fock 
approximationi^ for the intraband interaction (V^"^^). 
A more rigorous analysis requires solving the Bethe- 
Salpeter equation.— Moreover, we have normal-ordered 
the density operators (: '^l^c ■= *l*c and : ^I'^v := 
-•^v'i'l), and have used *c|0) = = *J,|0). 

The interband Coulomb matrix element in Eq. ([6]) 
reads 
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F(k,k') 



2e|k-k' 



EkEk' 



1 + k • k' - iA (— + ^ ) £ • (k X k') 

EkEk' J \Ek Ek' 



r 



(7) 



where e is the momentum-dependent static dielectric con- 
stant at the surface of the TI. A proper microscopic the- 
ory of e requires the evaluation of a Lindhard functionii 
that includes contributions from both surface and bulk 
electronic bands. In lieu of a detailed knowledge of the 
full band structure, we hereafter regard e as an empir- 
ical, momentum-independent quantity. The assumption 
of momentum-independence will be justified a posteriori 
(Figs. ([2]) and ([3])), when we find that x(k) is peaked at 
/c ~ 0. Due to this, the most important Coulomb ma- 
trix elements in Eq. (jB]) carry a small momentum transfer 
(|k — k'l ~ 0) and accordingly the long- wavelength limit 
of the static dielectric constant furnishes a satisfactory 
description of screening. 

In Eq. ([7]), the second and last terms within the brack- 
ets are unconventional in the sense that they depend on 
the relative angle between the momenta of the interact- 
ing particles. In particular the last term is chiral, i.e. 
sensitive to the relative angular momentum of the inter- 
acting particles, and arises purely due to the broken time- 
reversal symmetry (it changes sign under A — >■ — A.) For 
A = the Coulomb matrix element agrees with that of 
graphene. Note also that V(k,k') =^ V*{k',k). Eq. ^ is 
the two-particle analog of the skew scattering that con- 
tributes to the anomalous Hall effect in ferromagnetic 
metalsii^ In that case the probability of scattering off 
a nonmagnetic impurity is dependent on the chirality or 
handedness between the initial and final momenta. In our 
case, the strength with which an electron and a hole in- 
teract is skewed towards a particular sign of their relative 
angular momentum. This has important implications for 
the exciton binding energies, as we discuss below. 

Continuing with Eq. ([5]), the renormalization of the 
band-gap in the Hartree-Fock approximationi^ is given 
by the difference between self-energies of the conduction 
and valence band: 



£2 

i< 

= -y 



^ elk- 

k' ' 



-(|(c,kKk')p-|(^.,kKk')P) 



27re2 A2 ^ h^v^k ■ k! 



e|k-k' 



EkEk' 



(8) 



where a momentum cutoff kc must be introduced in 
order to regularize the ultraviolet divergence that re- 
flects the breakdown of the low-energy effective the- 
ory. (5Sk does not depend on the direction of k, and 
thus preserves the rotational symmetry of the problem. 
i5Sk=o = 2gA\n(2kcHvp / A) renormalizes the gap at the 
Dirac point, and as expected it vanishes for A = 0. For 
k ^ A/{hvp), SJ^k is generally nonzero even if A — >■ 0. 
In order to capture the physical role of the self energy, it 



is illustrative to represent it in operator language, 

k 

-Eii(i,.<ij-('^^+"--'')(i: 

k 



which indicates that Coulomb interactions increase the 
bare Fermi velocity and the bare magnetic gap via 



VF vpil + 5^k/2Ek) 
A ^ A{1 + S^k/2Ek). 



(9) 



For k ^ A/{hvp) this enhancement amounts to up to 
10% of the bare value. Eq. ^ suggests an explanation for 
why the Fermi velocity of the topological surface states 
in Bi2Se3 has been found to decrease^^ under Cu doping. 
In effect, the Coulomb interaction becomes weaker due 
to increased e and onset of screening by doped carriers, 
and therefore the renormalized Fermi velocity is smaller 
than for the undoped case. 



A. Bound Electron-Hole Pairs 

Eq. ([5]) is a hydrogen-like Schrodingcr equation for an 
electron-hole pair with zero COM momentum. As such 
we expect a solution that will consist of a discrete set of 
bound states (excitons), as well as a continuum of scat- 
tering states. Because of the unconventional momentum 
dependence of the Coulomb matrix elements in Eq. ([7|), 
Eq. ^ must be solved numerically. We begin by decom- 
posing the wavefunction into angular momentum chan- 
nels x(k) = Xm(fc) cxp(im0), which leads to 



i2Ek+SJ:k)Xm{k)- 



°° dk'k' 



where 



27r 



,fc')-^'^ge""V(k,k') 



Vm{k,k')Xra{k') = ephXm(fc), 

(10) 



V„^{k 



(11) 



and we have used (l/L^) Ek ^ / d^/i'^'^f- In Eq. ^ 
= 0k — 0k' is the angle between k and k'. While 
Eq. (|lip is physically meaningful only for k,k' < k^ the 
rapid decay of the electron-hole wavefunction in momen- 
tum space (see below) justifies taking the entire range 
of momenta. It is clear from Eq. ([TU|) that m is a good 
quantum number for electron-hole pairs with vanishing 
COM momentum. This is the U(l) symmetry identified 
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in the previous section, and accordingly m is the projec- 
tion of the total angular momentum along the z-axis . 
We solve Eq. ([T0| numerically using a modified quadra- 
ture method^S that carefully treats the singularity of the 
Coulomb matrix element at k = k' (see Appendix). 



3p 
2.5- 
2- 



renormalized bandgap 
»m=-l 
im=0 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



FIG. 1; Excitation energies as a function of the dimension- 
less Coulomb coupling g = {e^ / e) / {2hvF), for an electron-hole 
pair located at the surface of a TI with a magnetic gap. m 
is the quantum number for the projection of the total angu- 
lar momentum along the z-axis and Eg = 2A -I- (5Ek=o is the 
renormalized energy gap. As derived in the text, the renor- 
malized gap grows linearly with g, and its precise magnitude 
depends logarithmically on an ultraviolet momentum cutoff. 
The exciton binding energy is given by Eg — £ph and is found 
to obey Eq. (|13|) . For each m, only the lowest-lying bound 
state is shown. In the weak couphng regime (which is rele- 
vant for most TIs discovered to date) the m = — 1 exciton 
has the largest binding energy. However, the m = channel 
is most prone to an excitonic instability (at g ~ 0.4). The 
critical value of g for the instability can be roughly estimated 
from Eq. (|13|) via eb{g — >■ gc) — 2A, or else from dimensional 
arguments outlined in the text. 




FIG. 2; Momentum space wave function for the lowest-lying 
electron-hole pair with m = 0. as = hvF/{Ag) is an effec- 
tive Bohr radius (the effective mass in the paraboUic band 
approximation is A/2vp). The wavefunction (normalized 
through Jg°° d«:«:|xm(K)P = 1, where k. = kas) is peaked at 
kas ~ 0{i) {hvpk/A ~ (J <C 1), and its spread in momentum 
space increases as interactions become stronger. Inset: the 
z-component of the spin density for the same exciton state. 
Weakly interacting particle-hole pairs are spin polarized along 
the direction normal to the surface. This spin projection de- 
creases abruptly in the vicinity of the excitonic instability. 



excitons with a given value of the angular momentum to: 



Eo 



(|to| + i/2y 



(ordinary 2D excitons). (12) 



The decrease in the binding energy of the strongest 
bound exciton with m may be interpreted by invoking 
a centrifugal potential that reduces the effective attrac- 
tion between the electron and the hole for higher angular 
momentum channels. In contrast to Eq. p2p . our numer- 
ical analysis for excitons on topological surfaces gives 



En 



In Fig. ([T|) we plot the energy of the lowest lying ex- 
citon states for different values of m. At weak Coulomb 
coupling the binding energy is largest for m = — 1, and 
the excitation energies for m = 1 and to = — 1 are 
different. Such asymmetry in the sign of the angular 
momentum stems from the aforementioned chirality of 
the Coulomb matrix clement (Eq. (O), which enhances 
the electron-hole attraction for angular momenta aligned 
with the Zccman field Az. As such, to and —to arc in- 
terchangeable via A — > — A. This is all in stark contrast 
to conventional hydrogen-like excitons. In ordinary two- 
dimensional excitons with a reduced effective mass TOeff, 
the binding energy is given by Cb = Eo/{n + l/2)'^, where 
n = 0,1,2,..., Eq = (TOeff/TOee^)13.6eV and TOe is the 
bare electron mass. It is instructive to write down the 
expression for the largest binding energy of conventional 



b.m 



(|to + 1| + 1/2)2 



(topological 2D excitons), 



(13) 

where Eq = aAg^ and a is a constant of order one. Ag^ 
is the Coulomb energy scale corresponding to electron- 
hole pairs with a reduced effective mass A/2vp. This 
mass agrees with the inverse of the non-interacting band 
curvature at k = 0. The remarkable difference between 
Eq. (|12p and Eq. (|13p is a shift in the angular momen- 
tum quantum number (i.e. to — >■ to -I- 1), which is con- 
sistent with the nontrivial Berry phase of the topological 
surface states. Given their maximal binding energy for 
to = —1 as well as their asymmetry under to — >■ —to, 
the excitons occurring at topological surfaces with mag- 
netic order can be described as "p-wave" and "chiral". 
As we show below, these unconventional features may be 
optically detectable. 
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FIG. 3: Momentum space wave function for the lowest-lying 
electron-hole pair with g — 0.1 and m = 1,0,-1. The m — 
— 1 exciton wavefunction is more spread out in momentum 
space, which is partly why the corresponding binding energy 
is larger. 



For g ~ 0.1 and A ~ 20nieV, the binding energy for 
the lowest-lying exciton is about ImcV. Though small, 
this binding energy is comparable to that of biexcitons in 
GaAs quantum wells, which have been seen in photolumi- 
nescence experiments.— In the next subsection we show 
that excitonic features in the optical absorption spectrum 
can be noticeable even for g ~ 0.05 and A ~ lOmeV. 

For larger Coulomb coupling (g > 0.4), the binding en- 
ergy for the m = exciton surpasses that of the m = — 1 
exciton, and leads to an excitonic instability towards an 
"s-wave" Bose condensate. The critical value of g for the 
onset of the excitonic instability can be estimated from 
dimensional analysis. At zero temperature and when 
the Fermi energy is in the gap, the only lengthscale^^ 
for non-interacting electrons is = hvp/^, whereas 
the only lengthscale related to Coulomb interactions is 
= e^/(eA). The validity of our calculations is lim- 
ited to the weakly interacting regime where ^ ■^A- In 
contrast, when 3> the single-exciton approxima- 
tion breaks down and the physics is governed by many- 
body interactions. Therefore the excitonic instability oc- 
curs at — Ca, which yields gc — 0.5. This value 
agrees well with the numerical results of Fig. ([T]). For 
Wi? ~ 4 X lO^m/s, g = 2.72/e and thus the excitonic in- 
stability is reached at e ~ 7. At present we are not aware 
of any TI that has such a small value of the dielectric 
constant at its surface. 

Fig. illustrates the wavefunction for the m = 
exciton in momentum space. Most of its weight is con- 
centrated at kuB ^ where as = hvp/(Ag) is an 
effective Bohr radius. For strongly interacting electron- 
hole pairs, the wavefunctions become more delocalized 
in momentum space. This trend reflects the tendency of 
the electron and the hole to stay close to one another 



in real space, so as to optimize the Coulomb interac- 
tion. The inset of Fig. ^ shows the spin density of the 
TO = exciton. In the non-interacting case, the spinor 
of the conduction band at momentum k is perfectly an- 
tialigned with respect to that of the valence band. There- 
fore, non-interacting particle-hole excitations are spinful. 
The X- and y-components of the spin density vanish by 
symmetry, while the z-componcnt of the spin density for 
electron-hole pair with angular momentum m is 

1 r dkk A 

= -(xc|a|+al|xc) = ^ / __|x„„(fc)p, (14) 

where n labels the different eigenstates for a given value 
of TO. As evidenced by Fig. Mz varies little as a 
function of g up until g g,.. It is sensible that the spin- 
polarization of electron-hole pairs be largely immune to 
weak Coulomb interactions, because the latter are spin- 
rotationally invariant. 

In Fig. ([3]) we plot the momentum space wave func- 
tions for the TO = 1,0,-1 excitons. The to = — 1 
(to = 1) eigenstate is most (least) delocalized in mo- 
mentum space. This is a consequence of the chirality in 
the Coulomb interaction, which makes the electron-hole 
attraction largest for to ~ — 1 and smallest for to — 1. 
Hence the to = — 1 wavefunction is more delocalized in 
momentum space. 

B. Optical Absorption 

The knowledge of particle-hole excitation energies and 
wavefunctions at P = provides a gateway for deter- 
mining the optical conductivitj«2^i2i of topological surface 
states whose chemical potential is placed inside a mag- 
netic energy gap. From Fermi's golden rule, the optical 
absorption for photons of frequency u is 

= Y ^ '^"'^^ ■ ^''^^'''^(^h - M, (15) 

a 

where a is the set of quantum numbers describing a 
particle-hole excitation and |0) is the ground state in 
absence of particle-hole excitations. Ignoring electric 
fields (which are included in the next section) we have 
a — {71,111), where to refers to the angular momentum 
quantum number and n counts the different eigenstates 
for a given to. v = dHo/hdlL ~ vp{(T x z) is the ve- 
locity operator for non-interacting electrons (we neglect 
the rcnormalization of the Fermi velocity due to interac- 
tions), A is the vector potential and 

(«|i;^|0) = ^xr™(fc)e™*(ck|i;>k) (16) 

k 

is the electric dipole transition matrix element. 

Because [Vj'Hq] 0> the velocity operator induces in- 
terband transitions and is able to produce magnetized 
excitons. The optical absorption vanishes for A = Az 
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and is symmetric in the xy plane. For linearly polarized 
light (e.g. A = Ax, though one can immediately ver- 
ify that any linear polarization in the xy plane yields an 
identical optical absorption), the dipole matrix element 
is 



(ck|z;^|wk) = 



'IVF 



e*^'' cos^ 



2 



sm 



2 t'k 



(17) 

where 6'k = cos ^(A/iJk) and (/3k = (/)+7r/2 are the spinor 
angles introduced in Section III. It follows from Eqs. (fT6|) 
and (jl7p that only m = 1 and m = ~~1 states arc optically 
active for linearly polarized light, whereas the m = 
exciton is dark. This selection rule differs from ordinary 
hydrogen-like excitons and even from excitons in many 
topologically trivial spin-orbit coupled systems, wherein 
m = is a bright exciton. It is likewise different from 
the optical selection rules in gated bilayer graphenej^ 
For circularly polarized light (A = A{x ± iy)/^/2) we 
have 



(ck|w+|wk) 
(ck|?;~ |uk) 



V2wj?e''^'' cos^ 



-\/2u'_Fe sm 



2 

2^ 

2 ' 



(18) 



where v'^ = {v^ ± iw^)/-\/2. Thus only the m = —1 
(m = 1) exciton is optically active for right- (left-) cir- 
cularly polarized light. In Figs. ^ and ([5|) we plot the 
optical absorption for linearly and circularly polarized 
light. We incorporate the finite lifetime of the excitons 
by broadening the Dirac delta function 5{en,m — fvuj) into 
a Lorcntzian 



hjj) 



7/2 



(n,m) 

ph 



/iw)2 +72/4 



(19) 



where 7 is the linewidth. We have taken 7 ~ 0.6mcV, 
which is fairly small but commonly obtained in a variety 
of semiconducting quantum wellsj2£ 

In the non-interacting case the optical absorption in- 
creases linearly at frequencies above the gap; this is sim- 
ply a manifestation of the Dirac-like density-of-states. 
For frequencies in the vicinity of the gap, electron-hole in- 
teractions enhance the optical absorption. When the sur- 
face is irradiated with either linearly polarized or right- 
circularly polarized light, an excitonic absorption peak is 
visible even in samples with modest gaps (A ~ lOmcV) 
and weak interactions [g ~ 0.05 or e ~ 50). Part of the 
reason for this prominent excitonic peak is that the "co- 
herence factor" cos^(0/2) in Eqs. pT)) and ([T^ has most 
of its weight at hvpk <C A, which is precisely the region 
of momentum space where the exciton wavefunction is 
peaked (recall Fig. 

In case of left-circularly polarized light, the optically 
bright m = 1 excitons are esentially merged with the 
continuum at experimentally relevant Coulomb interac- 
tion strengths (recall Fig. ([T|)), and moreover the coher- 
ence factor sin^((?/2) in the second line of Eq. is 
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FIG. 4: Optical absorption for linearly polarized light, as 
a function of the dimensionless Coulomb coupling A = 
20meV, 7 = O.GmeV, and Ro is the absorption coefficient for 
free electron-hole pairs aX u = 2A. Eg = 2A -|- (5Ek=o is 
the renormalized energy gap. Most of the subgap absorption 
originates from m = — 1 (or m = 1 if one reverses the direc- 
tion of magnetization). Inset: Optical absorption spectrum 
for linearly polarized light, with g — 0.1, A = lOmeV and 
7 = 1.2meV. The intensity of the m = —1 intensity peak is 
reduced but still noticeable. 



largely depleted in the region of momentum space where 
the excitonic wavefunction is localized. Consequently, 
the absorption spectrum of left-circularly polarized light 
shows no significant excitonic peaks and the absolute in- 
tensity of absorption near the bandgap is much smaller 
(Fig. (O). Such a large asymmetry between left- and 
right-circularly polarized light should be experimentally 
testable. Keeping in mind that right- and left-circularly 
polarized absorption spectra map onto one another un- 
der A — > — A, one simple way to detect the asymmetry 
between them is by shining circularly polarized light and 
then switching the direction of magnetization at the sur- 
face of the topological insulator by 180°. 

We finalize this section with a few caveats on how to fa- 
cilitate the experimental detection of the aforementioned 
phenomena. On one hand, the magnetic overlay must be 
a charge insulator with an optical gap that is large com- 
pared to the energy gap of the surface states. This re- 
quirement is clearly fulfilled in magnetic ions (e.g. Fe"^"*", 
Mn^+) or in magnetic semiconductors such as21 EuO, 
which have O(leV) bandgaps. On the other hand, light 
can excite spatially uniform magnons in the ferromag- 
netic overlay. However, the typical energy scale for these 
collective cxcitations^Si is 0.05meV ^ 2A, and therefore 
the magnon absorption lines will be well separated from 
excitonic features. Perhaps the most convenient system 
for testing our predictions is Bi2Se3 with Fe and Mn ions 
doped in the bulk or at the surfacei^ Bulk Mn doping 
leads to A '-^ 7meV at 1% doping; substantially larger 
magnetic gaps have been reported through increase in 
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FIG. 5: Optical absorption for circularly polarized light, 
g = 0.12, A = 15meV and 7 = 0.6meV. i?o is the ab- 
sorption coefficient for <; = and uj = 2A, taken as the 
average of left-circularly and right-circularly polarized lights. 
Eg = 2A + 5Sk=o is the renormalized energy gap. For right- 
circularly polarized light the most strongly bound exciton 
(m = —1) is optically bright. For left-circularly polarized 
light m = —1 is optically dark and the absorption spec- 
trum is missing a prominent excitonic peak. The right- and 
left-circularly polarized spectra map onto one another under 
A — ^ —A, namely under switching of the magnetization. In- 
set: Zoomed in absorption spectrum for the left-circularly 
polarized light. The coherence factor sin^(S/2), mentioned in 
the text, diminishes the absorption intensity considerably. 



the dopant concentration. 



IV. ELECTRON-HOLE PAIRS ON A 
ONE-DIMENSIONAL SURFACE 



where fc is the momentum along the axis of the cylin- 
der, m is the azimuthal projection of the total angular 
momentum, v is the magnetic flux (in units of the flux 
quantum) that threads the cross section of the nanowire, 
and ro is the radius of the wire. For = 0, the en- 
ergy spectrum at the surface is gapped and all energy 
levels are doubly degenerate (Ek,m=-\ = Ek,m=o)- For 
V = 1/2, m = constitutes a gapless, non-degenerate, 
topological surface state. When 1/ & (0, 1/2) the surface 
spectrum is gapped but the lowest-lying energy states 
are non-degenerate; this is the situation for which the 
calculation of the present section applies. It follows from 
Eq. (|20p that the bandgap for non-interacting electrons 
is 2 A = {2hvF /ro)]!^ — 1/2 1, which is tunable by an ex- 
ternal magnetic field along the axis of the cylinder. An- 
other way to open a uniform energy gap is to coat the 
nanowire with a radially magnetized ferromagnet, though 
this is experimentally cumbersome, i/ ~ 1/2 is achievable 
with magnetic fields B < lOT even for radii as small as 
To ~ lOnm. When ly ~ 1/2, it is a good approximation 
to focus on the to = state only, because it is well sep- 
arated (at least near k = 0) from higher energy states. 
For instance, Ek=o,m=-i = 3Ek=o.,m=o for = 1/4. For 

~ 1/4 and rg ~ lOnm, A ~ O(lOmeV) is comparable 
to the magnitude of the exchange field considered in the 
previous section. 

When £ ^ the problem at hand contains no sym- 
metries other than the conservation of the COM canon- 
ical momentum P. Hence it is more convenient to solve 
the two-body Dirac equation in real space. Our starting 
model for the surface of a TI nanowire with radius tq and 
length L is 

n = -VF {al - a^) P+^{cT, + an) ■ {-px + —(l) 

2 \ VF J 

- -jr^ r-efy, (21) 



Thus far we have considered two dimensional surfaces 
with translational and rotational symmetry. In this sec- 
tion we concentrate on a one-dimensional variant of the 
interacting two-body problem (Eq. ([T|)), which is perti- 
nent to the surface of a quasi-one-dimensional TI. The 
motivation for this investigation lies in the well-known 
fact that excitonic effects are stronger in reduced dimen- 
sions. Besides, the one-dimensional model is numerically 
simple enough to allow us to learn with relative ease how 
excitonic effects are altered under external electric fields 
as well as in presence of a finite COM momentum for the 
electron- hole pairs. In regards to the latter, some of the 
results derived in this section (e.g. Eq. ([221) ') ^.re relevant 
for two dimensional surfaces. 

The low-energy band structure for the surface states 
of a cylindrical TI nanowire^ is 

E,,^^±nvF,l^^^^^±^^. (20) 



where we have defined the axis of the cylinder to be 
the y-axis and y e {~L,L) is the relative coordinate 
along the wire. The regularized form of the Coulomb 
interaction in Eq. (|2ip can be derived by averaging the 
three-dimensional Coulomb potential with radial enve- 
lope functionsi^ is a positive fitting parameter that 
takes the value of 0.3 in certain idealizations of wires. 

We solve Eq. (PT|) as an eigenvalue problem in real 
space, with the boundary conditions chosen so that the 
wavefunction vanishes at the ends of the wire. Starting 
from the non-interacting problem with P = £ = e^/e = 
0, we keep only a quarter of eigenvalues and eigenstates. 
These correspond to particle-hole pairs with positive ex- 
citation energies eph > 2 A, namely a hole in the valence 
band and an electron in the conduction band. There- 
after we analyze how these energy levels (and eigenstates) 
change under Coulomb interactions, electric fields and 
COM momenta. The eigenvalues and eigenstates we dis- 
card describe (i) a valence band hole and a valence band 
electron, (ii) a conduction band electron and a conduc- 
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tion band hole, 
band electron. 



fiii) a conduction band hole and a valence 




FIG. 6: Lowest-lying particle-hole excitation energy for a TI 
nanowire of radius ro = lOnm, length L = lOOnm and tivp ~ 
SOOmeVnm. In our numerical calculation, each discrete step 
along the wire equals 1 nm. The magnetic flux along the axis 
of the wire isv = 1/4 in units of the flux quantum, which leads 
to a non-interacting bandgap of 2A = 15meV. For smaller 
values of the fitting parameter 77, the exciton binding energies 
become larger. 



In Fig. ([6]) we plot the lowest lying particle-hole excita- 
tion energy as a function of g, for £ = ilx = = P = 0- 
Narrower wires show larger exciton binding energies: for 
wires of radius ~ lOnm they may be ~ 3 — 5 times larger 
than the binding energies of the two dimensional excitons 
studied in the previous section. 

In Fig. ([7]) we display the influence of P 7^ and radial 
exchange fields in the particle-hole excitation energy (az- 
imuthal exchange fields amount to a shift in p and arc in- 
consequential). The binding energy has a non-monotonic 
dependence on P; this is markedly different from con- 
ventional systems wherein the COM and relative coordi- 
nates are decoupled and thus the exciton binding energy 
is independent of P (at least in absence of quantizing 
magnetic fields). For P <C ^/vp, excitons are formed 
predominantly out of electrons (holes) located near the 
bottom (top) of the non-interacting conduction (valence) 
band. Accordingly we find eph(P) ~ eph(O) + P'^/{2M), 
where the exciton mass is given by AI = A/{2vp). As 
P increases, the particles and holes forming the exciton 
necessarily explore the linear part of the non-interacting 
energy-momentum dispersion. Thereby excitons become 
heavier; in fact for P ^ A/vp we obtain deph/dP ~ 2vf, 
which translates to infinitely massive particle-hole pairs. 

The observations from the preceding paragraph are 
also applicable to the two-dimensional magneto-excitons 
discussed in Section II. Yet unlike in nanowires, in two- 
dimensional topological surfaces the magnetic gap is gen- 
erated by exchange interactions. These exchange in- 




FIG. 7: Lowest-lying particle-hole excitation energy as a func- 
tion of the canonical COM momentum P. L — lOOnm, 
A_ = and g = 0.12. The dependence of the particle- hole en- 
ergy on P is quadratic at P ^ A/vp (light excitons) and be- 
comes linear as P 3> A/vp (heavy excitons). When A_ 7^ 0, 
the dispersion is shifted. A nearly identical dispersion is ob- 
tained (not shown) as a function of a radial exchange field hy 
that couples to + a^. The effect of this exchange field is 
to renormalize the bandap via A 



+ hy. Azimuthal 



exchange fields (which couple to +a^) amount to a redefi- 
nition of the origin in momentum space and do not change the 
particle- hole excitation spectrum. Inset: The exciton binding 
energy has a non-trivial dependence on P. 



teractions are generally asymmetric between electrons 
and holes. An interesting consequence of this asymme- 
try is that it introduces a mechanical COM momentum 
n = P — /vpzxn, where A_ denotes the difference in 
the exchange fields between electrons and holes (see Sec- 
tion II for the definition of A±). When A_ ^ 0, Fig. © 
remains valid insofar as P is replaced by 11 in the hori- 
zontal axis. As a result optically generated excitons (i.e. 
P = excitons) carry a nonzero COM group velocity in 
presence of in-plane exchange fields: 



lim 

P->0 



geph(P) 
dP 



Vp 



-z X VL. 



(22) 



Eq. ([^ was derived for A_|z x 0| ^ f2zA+, and yields 
Vg « vp- This nonzero group velocity of optically ex- 
cited excitons, combined with their nonzero spin polar- 
ization along the z-axis (see Eq. (jl4p and the discussion 
therein), opens prospects for creating cxcitonic spin cur- 
rents and spin Hall effects. Similarly, by tilting 17 away 
from z one can map the excitonic dispersion relation 
(Fig. 0) optically,^ 

Finally, we investigate the photoluminescence 
(Eq. (Uni) of the nanowire in presence of electric fields. 
Let us denote the eigenstates of Eq. (PT|) as x'^ (y)j 



where the index a is such that e. 



(a) 
ph 



> 2A in the 



non-interacting case (recall the discussion above), and 
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FIG. 8: Optical absorption (in arbitrary units) for non- 
interacting electrons illuminated by linearly polarized light 
(with its polarization vector along the axis of the nanowire). 
L — 150nm is the length of the nanowire, 2A — 15meV is 
the non-interacting bandgap, v — 1/4 is the magnetic flux 
threading its cross section, 77 — 0.3 is the fitting parame- 
ter that enters in the regularized Coulomb interaction, and 
7 — 1.2meV is the excitonic linewidth. In absence of electric 
fields, the absorption spectrum reflects the one-dimensional 
density of states (modulated by the "coherence factors" of 
Eq. (|23p ') at frequencies above the bandgap. In presence of 
electric fields the subgap absorption is enhanced and Franz- 
Keldysh oscillations appear. 



a (cr') is the spin label of the electron (hole). Then the 
dipole matrix element is given by 

(a|z;^|0)=^ ^ x7'ik)M.,>{k){ck\vy\vk), 

k aa'£{t,i} 

(23) 

where is the velocity operator along the wire, Xa{k) = 

J^^dye-^'^yXaiy) and Maa' = *e.(fc) ® ^ha'{-k) = 

(cos(0fc/2), - sm{0k/2)) ® (cos(0fc/2), - sm{0k/2)). The 
optical absorption is nonzero provided that the polar- 
ization of the impinging light has a nonzero projection 
along the axis of the nanowire. Fig. (|8]) shows that elec- 
tric fields enhance the subgap optical absorption due to 
particle-hole excitations with eph < 2A being formed be- 
tween different points along the wire. The electric-field- 
induced oscillations in photoluminescence, i.e. Franz- 
Kcldysh oscillationsj^i^ arc visible for electric fields 
£ > 250 V/cm. Fig. ([5]) illustrates that the influence of 
a given electric field in the excitonic absorption peak de- 
creases as the Coulomb interactions are made stronger. 
This is reasonable because the electron and the hole tend 
to stay close to one another in order to maximize their 
Coulomb attraction, whereas the electric-field term in the 
Hamiltonian vanishes when the electron and the hole are 
located at the same point in space. Incidentally, we find 
that an electric field induces no net spin polarization of 
the exciton. 



FIG. 9: Optical absorption for interacting electrons illumi- 
nated by linearly polarized light. The parameters of the 
nanowire are identical to those of the previous figure, with 
7 = 1 meV. Eg — 2A is the bandgap for non-interacting elec- 
trons. When £ = 500V/cm, the influence of the electric field 
in the optical absorption is clearly visible for g = 0.03 and 
less discernible for g — 0.12. 



V. CONCLUSIONS 

Motivated in part by the objective to understand the 
interplay between Coulomb and spin-orbit interactions 
on topological insulators with broken time reversal sym- 
metry, we have studied the excitonic photoluminescence 
on a magnetized surface of a strong TI. Throughout this 
study we have encountered a number of potentially in- 
teresting properties. First, optically excited electron- 
hole pairs are magnetized in the direction perpendicu- 
lar to the topological surface. Second, the relative mo- 
tion of an electron-hole pair is entangled to the ccnter-of- 
mass motion. Accordingly, the excitonic Berry phase is 
nonzero and a spin-Hall effect of excitons ensues. Third, 
Coulomb interactions enhance both the Fermi velocity 
and the mass of Dirac fermions. Fourth, in two dimen- 
sional topological surfaces, the combination of spin-orbit 
interaction and broken time reversal symmetry leads to 
a skewness or chirality in the effective Coulomb inter- 
action between an electron and a hole. The underly- 
ing cause for this chirality resides in the nonzero Berry 
phase of the topological surface states. Consequently, the 
strongest bound state occurs for an electron-hole pair 
whose total angular momentum has a projection h (or 
~h, depending on the direction of magnetization at the 
surface) along the normal to the surface. This "p-wave" 
and "chiral" exciton becomes optically bright upon shin- 
ing circularly polarized light with the proper handedness, 
and is optically dark upon shining of circularly polarized 
light with the opposite handedness. For fixed handedness 
of the impinging light, bright-to-dark transitions can be 
induced by switching the direction of magnetization by 
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180°. Fifth, so long as there is an electron- hole asymme- 
try in the exchange coupling between topological surface 
states and magnetic ions, optically excited electron-hole 
pairs have a nonzero center-of-mass group velocity that 
may be tuned by tilting the magnetization away from the 
normal of the surface. 

Many of the anticipated peculiarities may be ex- 
perimentally accessible in topological surfaces that 
host weak-to-moderate Coulomb interactions {g — 
{e'^/e)/{2hvp) > 0.1), promote sufficiently long exciton 
lifetimes, and have their chemical potential in the energy 
gap (both in bulk and at the surface). If the chemi- 
cal potential intersects the topological surface bands, an 
optical gap will certainly be still present. Nonetheless, 
excitonic effects will be substantially weakened because 
exciton binding energies^ and lifetimes^^ both decrease 
rapidly in the metallic regime. 

Aside from the task of characterizing magneto-excitons 
that form between Landau levels under perpendicular 
magnetic fields, we can envision various directions along 
which the work presented here can be expanded. On 
one hand, it will be interesting to consider ultrathin TI 
ribbons, wherein two Dirac cones (one on each surface) 
and spatially indirect excitons might contribute to pho- 
toluminescence. On the other hand, full-fledged elec- 



tronic structure calculations of the topological surface 
states will be crucial^ in order to make quantitative con- 
tact with optical measurements that are likely to be per- 
formed in the foreseeable future. 
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Note added in proof. - After submitting this manuscript 
we became aware of a recent preprint by A. Cortijo,— 
wherein an interesting chirality is unearthed in the pho- 
ton propagator on magnetized topological surfaces. The 
chiral terms in the photon propagator are relativistic and 
thus negligible in the context of our work, where we 
employ the dominant non-relativistic piece of the pho- 
ton propagator. The chirality of the interband Coulomb 
matrix elements in our Eq. ([7]) stems from the spin- 
momentum locking of the eigenstates. 
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Appendix A: Numerical Solution of Eq. (|10|) 



We begin by defining the dimensionless quantities k = hvpk/A, Vm{k,k') — Vm{k, k') / (Airhvp), = eph/(2A) 
and (5E = i5E/(2A), upon which Eq. pO)) can be rewritten as 

/"OO 

[(P + 1)1/2 ^ sf:]xmCk) - / dk'k'V^Ck, ~k')Xrn{k') = e~phx™(^). 



Next, we make a change of variables via k = tan(7ra;/2) and discretize x G [0, 1] in N segments, where the spacing 
between subsequent Xi and the corresponding weights Wi are chosen according to a A^-point Gauss-Legendre quadrature 
method. This yields^S 



^ Uk 



[{kj + 1)1/2 ^ g^^^^^ _ yU^^ _ ^ _ ^^^^^ ^ g^^^^ (Al) 



where = V;„(fci, kj) for j, 

r°° ^ ( dk\ 

Vj:^J^ dk'k'fCk,k')-Y. f{k,kykA^j w,, (A2) 



and /(fc. A:') = 2V^(fc, k')k'^/{k'^ + k 2). The first term in Eq. (|A2[) is evaluated numerically after shifting k' so as to 
remove the Coulomb singularity. In the numerical evaluation of Eq. (|A1I) we take N = 256; higher values of N do not 
produce any significant quantitative changes. 
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